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Abstract. We present a fully microscopic study of the ^^O+^^^Pb fusion using the density-constrained 
time-dependent Hartree-Fock theory. The calculated fusion cross-sections are in good agreement with the 
experimental data for the entire energy range indicating that the incorporation of dynamical effects is 
crucial in describing heavy-ion fusion. 
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1 Introduction 

Heavy-ion fusion reactions are a sensitive probe of the 
size, shape, and structure of atomic nuclei as well as the 
collision dynamics. Fusion studies using neutron-rich nu- 
clei are becoming increasingly available. Examples include 
the experiments with heavy neutron-rich ^'^^Sn beams on 
^■^Ni [fl[2], and the fusion measurements for the ^'*Ni-}-^'*Ni 
system [3] . What characterized these experiments was the 
measurement of the fusion cross-sections at very low bom- 
barding energies. Recently, fusion cross-sections for the 
i6Q_,_208p|3 system were also measured at extreme sub- 
barrier energies [1] . Another experimental frontier for fu- 
sion is the synthesis of superheavy nuclei in cold and hot 
fusion reactions [HEIITIIH] . 

To date quantitative theoretical studies of fusion cross- 
sections have largely been limited to phenomenological 
methods, mainly the coupled-channels approach. Specif- 
ically, various effects such as neutron transfer [5], com- 
pression potential [9], modifications of the inner turning 
points via the approximate inclusion of neck formation 
[T^ , and complex potentials [13^ were introduced, in addi- 
tion to standard channel couplings, to address the recent 
data at extreme sub-barrier energies. This is surprising 
since the ^^0-|-^*'^Pb system should be an ideal system for 
coupled-channel calculations. A number of papers have 
been devoted exclusively to the ^^0-)-^"^Pb system T^, 
[To] . This underscores the fact that while phenomenological 
methods provide a useful and productive means for quan- 
tifying multitudinous fusion data it is desirable to make 
contact with the microscopic theories of nuclear structure 
and reactions. 

In the absence of a true many-body theory for cal- 
culating sub-barrier fusion all practical calculations are 
reduced to the calculation of a potential barrier between 
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the interacting nuclei and a subsequent calculation of tun- 
neling through the barrier. Many of the phenomenological 
methods use frozen densities for each nuclei and ignore 
all of the dynamical effects as well as the formation of a 
neck during the nuclear overlap. It has been demonstrated 
that for deep sub-barrier energies the inner part of the po- 
tential barrier plays a very important role [TTKT^. While 
the outer part of the barrier is largely determined by the 
early entrance channel properties of the collision, the in- 
ner part of the potential barrier is strongly sensitive to 
dynamical effects such as particle transfer and neck for- 
mation. Recently, an inversion method has been devised 
jlOj to obtain inter-nucleus potentials directly from the 
fusion data for the ^^O-t-^^^Pb system. In their conclusion 
the authors propose that dynamical effects, such as neck 
formation and coordinate dependent mass may be respon- 
sible for the modification of the ordinary potentials used 
in the coupled-channel calculations. 

It is generally acknowledged that the time-dependent 
Hartree-Fock (TDHF) theory provides a useful foundation 
for a fully microscopic many-body theory of low-energy 
heavy- ion collisions [14] . The success of the TDHF method 
is predicated on the expectation that the Pauli princi- 
ple plays an important role in simultaneously building up 
a time-dependent mean-field and suppressing the propa- 
gation of the strong N — N interaction terms. With the 
advent of computer technology it has now become fea- 
sible to perform TDHF calculations with no symmetry 
assumptions regarding the collision process and the effec- 
tive interaction, and with much more accurate numerical 
methods (TS]. At the same time the quality of effective 
interactions has also been substantially improved [IS]. In 
this sense, a description of the heavy-ion collisions via the 
TDHF method may provide a way to incorporate the dy- 
namical effects such as mass transfer and neck formation 
into the calculation of the ion-ion potential. 
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In Section 2 we outline the main features of one 
such method, the density-constrained TDHF (DC-TDHF) 
method |T7]. The DC-TDHF method was first used to 
study the fusion reactions of the ^^Ni-l-^'^^Sn system [HI 
[T^ . A re- measurement of the lowest energy data point [20] 
has moved this point exactly onto the DC-TDHF result. 
We have also performed calculations for the ^'*Ni-|-^'*Ni 
system i21, and showed that an excellent agreement with 
data is possible. As the DC-TDHF method contains no 
parameters or normalization its primary success depends 
on the quality of the description of the ground state prop- 
erties within the Hartree-Fock framework and the descrip- 
tion of the early stages of the ion-ion reaction by the 
TDHF theory. In this sense we expect the ^^0-|-^°®Pb sys- 
tem to be a good test since most Skyrme parametrizations 
include these nuclei in the fitting procedure. In Section 2 
we also discuss the calculation of ion-ion separation dis- 
tance, coordinate-dependent mass, and calculation of fu- 
sion cross-sections. In Section 3 the calculations for the 
i6Q_|_208p|-j syg^g]-|2 are presented. Final conclusions are 
given in Section 4. 

2 Formalism 

2.1 DC-TDHF method 

Recently, we have developed a new method to extract ion- 
ion interaction potentials directly from the TDHF time- 
evolution of the nuclear system. In the DC-TDHF ap- 
proach [TT] the TDHF time-evolution takes place with no 
restrictions. At certain times during the evolution the in- 
stantaneous density is used to perform a static Hartree- 
Fock minimization while holding the neutron and pro- 
ton densities constrained to be the corresponding in- 
stantaneous TDHF densities. In essence, this provides us 
with the TDHF dynamical path in relation to the multi- 
dimensional static energy surface of the combined nuclear 
system. The advantages of this method in comparison 
to other mean-field based microscopic methods such as 
the constrained Hartree-Fock (CHF) method are obvious. 
First, there is no need to introduce artificial constrain- 
ing operators which assume that the collective motion 
is confined to the constrained phase space, second the 
static adiabatic approximation is replaced by the dynam- 
ical analogue where the most energetically favorable state 
is obtained by including sudden rearrangements and the 
dynamical system does not have to move along the val- 
ley of the potential energy surface. In short we have a 
self-organizing system which selects its evolutionary path 
by itself following the microscopic dynamics. Some of the 
effects naturally included in the DC-TDHF calculations 
are: neck formation, mass exchange, internal excitations, 
deformation effects to all order, as well as the effect of nu- 
clear alignment for deformed systems. In the DC-TDHF 
method the ion-ion interaction potential is given by 

V{R) = Edc{R) ~ Ea, ~ Ea, , (1) 

where i?DC is the density-constrained energy at the instan- 
taneous separation R{t), while and are the bind- 



ing energies of the two nuclei obtained with the same effec- 
tive interaction. In writing Eq. ([!]) we have introduced the 
concept of an adiabatic reference state for a given TDHF 
state. The difference between these two energies repre- 
sents the internal energy. The adiabatic reference state is 
the one obtained via the density constraint calculation, 
which is the Slater determinant with lowest energy for the 
given density with vanishing current and approximates the 
collective potential energy [22]. We would like to empha- 
size again that this procedure does not affect the TDHF 
time-evolution and contains no free parameters or normal- 
ization. 

In addition to the ion-ion potential it is also possi- 
ble to obtain coordinate dependent mass parameters. One 
can compute the "effective mass" M{R) using the conser- 
vation of energy 

where the collective velocity R is directly obtained from 
the TDHF evolution and the potential V{R) from the den- 
sity constraint calculations. 

2.2 Calculation of R 

In practice, TDHF runs are initialized with energies above 
the Coulomb barrier at some large but finite separation. 
The two ions are boosted with velocities obtained by as- 
suming that the two nuclei arrive at this initial separation 
on a Coulomb trajectory. Initially the nuclei are placed 
such that the point x = in the x — z plane is the center 
of mass. During the TDHF dynamics the ion-ion separa- 
tion distance is obtained by constructing a dividing plane 
between the two centers and calculating the center of the 
densities on the left and right halves of this dividing plane. 
The coordinate R is the difference between the two cen- 
ters. The dividing plane is determined by finding the point 
at which the tails of the two densities intersect each other 
along the x-axis. Since the actual mesh used in the TDHF 
calculations is relatively coarse we use a cubic-spline inter- 
polation to interpolate the profile in the x-direction and 
search for a more precise intersection value. This proce- 
dure has been recently described in Ref. ^23J in great de- 
tail. We have confirmed with the authors of Ref. [23] that 
our R values are in exact agreement with their calculations 
for the i60-h^°^Pb system. 

The above procedure works well after the formation of 
the neck between the two ions. However, at some point 
when the two ions are substantially overlapping the pro- 
cedure starts to fail. While this does not occur for the 
R values used in this calculation, it may be desirable to 
find an alternate way of determining the R value. We 
have found that if one defines the ion-ion separation as 
R = Ra\/\Q2Q\, where Q20 is the mass quadrupole mo- 
ment for the entire system, calculated by using the colli- 
sion axis as the symmetry axis, and i?o is a scale factor 
determined to give the correct initial separation distance 
at the start of the calculations. Calculating R this way 
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yields almost identical results to the previous procedure 
until that procedure begins to fail and continues smoothly 
after that point. We have compared the two methods for 
the ^^O+^^O system for which the dividing plane is es- 
sentially at the origin for all R values. We find that the 
quadrupole method gives the correct R value even after 
substantial overlap. Of course the minimum value of R 
calculated this way is never zero but is determined by the 
quadrupole moment of the composite system. In this sense 
the above approach may be useful for fission studies. 

2.3 Fusion cross-section 



Integration of Eq. ^ yields 

R = fiR) ^ R = f-\R). (10) 

As a result of this point transformation, both the clas- 
sical Hamilton function, Eq. ([¥]), and the corresponding 
quantum mechanical Hamiltonian, Eq. ([5]), now assume 
the form 

H{R,P) = ^^+U{R) , (11) 

and the scaled heavy-ion potential is given by the expres- 
sion 

U{R)^ViR)^Vif-\R)) . (12) 



We now outline the calculation of the total fusion cross 
section using an arbitrary coordinate-dependent mass 
M{R). Starting from the classical Lagrange function 



1 



L(i?, R) = -M{R)R^ - V{R) , 
we obtain the corresponding Hamilton function 

p2 



H{R, P) 



2M{R) 



V{R) 



(3) 



(4) 



where the canonical momentum is given by P = M(R)R. 
Following the standard quantization procedure for the ki- 
netic energy in curvilinear coordinates [25] 



T = 



9 



._d_ 



-9^9' 



d 



(5) 



where g^u{q) denotes the metric tensor and g^'^(q) the 
reciprocal tensor, one obtains the quantized Hamiltonian 



H{R,P) 



M{R)-^PM{R)-^P +V{R) . (6) 



with the momentum operator P 
fusion cross cross-section 



'ihd/dR. The total 



(7) 



L=0 



can be obtained by calculating the potential barrier pen- 
etrabilities Tl from the Schrodinger equation for the rela- 
tive motion coordinate R using the Hamiltonian ^ with 
an additional centrifugal potential 



H{R,P) 



2M{R)R^ 



MR) = 



(8) 



Alternatively, instead of solving the Schrodinger equa- 
tion with coordinate dependent mass parameter M{R) 
for the heavy-ion potential V{R)^ we can instead use the 
constant reduced mass ji and transfer the coordinate- 
dependence of the mass to a scaled potential U{R) using 
the well known coordinate scale transformation HU 



3 Results 

In our mmierical calculations for the ^^O+^^'^Pb system 
we have chosen a Cartesian box which is 60 fm along the 
collision axis and 30 fm in the other two directions. The 
two nuclei are placed at an initial separation of 24 fm. 
Calculations are done in 3-D geometry and using the full 
Skyrme force (SLy4) [1^ as described in Ref. The nu- 
merical accuracy of the static binding energies and the de- 
viation from the point Coulomb energy in the initial state 
of the collision dynamics is on the order of 50 — 150 keV. 
We have performed density constraint calculations at ev- 
ery 10 fm/c interval or a total of 40 calculations. The 
accuracy of the density constraint calculations are com- 
mensurate with the accuracy of the static calculations. 




12_ 14 
R,R (fm) 

Fig. 1. (Colour on-line) Potential barriers, V{R), obtained 
from density constrained TDHF calculations using Eq. ([TJ at 
three different energies, Ec.m. ~ 76 MeV (black solid curve), 
E^.m. = 78 MeV (blue solid curve), and _Bc.m. = 90 MeV (red 
solid curve) . The three dashed curves correspond to the trans- 
formed potential in Eq. (|12p using the coordinate dependent 
masses. Also shown is the point Coulomb potential. 
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In Fig. [T] we show the ion-ion potentials obtained us- 
^ ' ing Eq. ^ calculated at three different energies, Ec.m. ~ 
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76 MeV (black solid curve), ^c.m. = 78 MeV (blue solid 
curve), and i^c.m. 90 MeV (red solid curve). The same 
energy dependence was also observed in the DD-TDHF 
calculations of Ref. [23] and seems to be more prevalent for 
heavier systems. For example in -'^^O+^^O system we saw 
a very small energy dependence even though the energy 
was tripled from 12 MeV to 36 MeV. As can be expected, 
at lower energies we observe more dynamical effects such 
as mass transfer, whereas the high energy results seem to 
approach the frozen-density limit [23' . This is also appar- 
ent in the coordinate dependent mass of Eq. Q shown 
in Fig. [2] for three different energies, Ec.m. — 76 MeV 
(black solid curve), Ec.m. = 78 MeV (blue solid curve), 
and Ec.m. = 90 MeV (red sohd curve). The i?-dependence 
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Fig. 2. (Colour on-line) Coordinate dependent mass M{R) 
scaled by the constant reduced mass fi, obtained from Eq. ([2]), 
at three different energies, Ec.m. = 76 MeV (black solid curve), 
-Ec.m. = 78 MeV (blue solid curve), and Ec.m. = 90 MeV (red 
solid curve). 

of this mass at lower energies is very similar to the one 
found in CHF calculations [24]. On the other hand, at 
higher energies the coordinate dependent mass essentially 
becomes flat, which is again a sign that most dynamical ef- 
fects are contained at lower energies. The peak at small R 
values is due to the fact that the center-of-mass energy is 
above the barrier and the denominator of Eq. ^ becomes 
small due to the slowdown of the ions. We have used the 
coordinate dependent masses shown in Fig. [2| to obtain 
the scaled potentials U{R) of Eq. ([T2|). These potentials 
are shown as the dashed curves in Fig. [TJ As we see the 
coordinate dependent mass only changes the inner parts 
of the barriers for all energies. Furthermore, the effect is 
largest for the lowest energy collision and diminishes as 
we increase the collision energy. 

We have obtained the fusion cross sections by numeri- 
cal integration of the Schrodinger equation using the con- 
stant reduced mass ^ and the scaled potential U {R) . The 
fusion barrier penetrabilities T^^Ecm..) are obtained by 



numerical integration of the two-body Schrodinger equa- 
tion using the incoming wave boundary condition (IWBC) 
method. IWBC assumes that once the minimum of the 
potential is reached fusion will occur. In practice, the 
Schrodinger equation is integrated from the potential min- 
imum, i?minj where only an incoming wave is assumed, to 
a large asymptotic distance, where it is matched to in- 
coming and outgoing Coulomb wavefunctions. The barrier 
penetration factor, Ti(i?c.m.) is the ratio of the incoming 
flux at i?niin to the incoming Coulomb flux at large dis- 
tance. Here, we implement the IWBC method exactly as 
it is formulated for the coupled-channel code CCFULL 
described in Ref. |26]. This gives us a consistent way for 
calculating cross-sections at above and below the barrier 
energies. The resulting cross-sections are shown in Fig. [S] 




E (MeV) 
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Fig. 3. (Colour on-line) Total fusion cross section as a function 
of Ec.m.- Three separate theoretical cross section calculations 
are shown, based on the energy-dependent DC-TDHF heavy- 
ion potentials V[R) and coordinate dependent masses M{R) at 
energies Ec.m. = 76, 78, 90 MeV, from which the scaled poten- 
tials U{R) were obtained. The experimental data (filled circles) 
are taken from Ref. [4]. 



We observe that all of the scaled barriers give a very 
good description of the fusion cross-section at higher en- 
ergies suggesting these cross-sections are primarily deter- 
mined by the barrier properties in the vicinity of the 
barrier peak. Naturally, the extreme sub-barrier cross- 
sections are influenced by what happens in the inner part 
of the barrier and here the dynamics and consequently the 
coordinate dependent mass becomes very important. As 
we observe in Fig. [3| we get a very good agreement with 
experiment for the barriers obtained at the lowest two col- 
lision energies, whereas the Ec.m. — 90 MeV curve under- 
estimates the cross-section at lower energies. Although not 
shown in Fig. [3] the cross-sections obtained using the the 
unsealed potentials V{R) and a constant reduced mass fi 
also agree well with the data for Ec.m. > 77 MeV but 
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either significantly over-estimate or under-estimate the 
cross-section at lower energies. 



4 Conclusions 

In summary, we have used the fully microscopic DC- 
TDHF method to study i^O-h^^^Pb fusion. The standard, 
parameter free approach of DC-TDHF yields potential 
barriers that can accurately reproduce the fusion cross- 
sections. The DC-TDHF approach has now been applied 
to study a number of systems with very promising results. 
Our results also agree with the conclusion of Ref. [TD] that 
for the proper description of fusion cross-sections at deep 
sub-barrier energies the dynamical effects such as neck 
formation and mass transfer must be included to modify 
the inner part of the potential barrier. Also, interesting 
are the implications of the energy dependent barriers as 
was also discussed extensively in Ref. [23]. These findings 
underscore the fact that additional modifications needed 
for phenomenological methods to explain the fusion cross- 
sections of the ^^0-|-^*'^Pb system may largely be due to 
the inadequacy of the approximations made in treating 
the nuclear dynamics. 

This work has been supported by the U.S. Depart- 
ment of Energy under grant No. DE-FG02-96ER40963 
with Vanderbilt University. 
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